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We propose a scenario which can explain large lepton asymmetry and small baryon asymmetry 
simultaneously. Large lepton asymmetry is generated through the Affleck-Dine (AD) mechanism 
and almost all the produced lepton numbers are absorbed into Q-balls (L-balls). If the lifetime 
of the L-balls is longer than the onset of electroweak phase transition but shorter than the epoch 
of big bang nucleosynthesis (BBN), the large lepton asymmetry in the L-balls is protected from 
sphaleron effects. On the other hand, small (negative) lepton numbers are evaporated from the 
L-balls due to thermal effects, which are converted into the observed small baryon asymmetry by 

CM . virtue of sphaleron effects. Large and positive lepton asymmetry of electron type is often requested 

from BBN. In our scenario by choosing an appropriate flat direction in the minimal supersymmetric 
standard model, we can produce positive lepton asymmetry of the electron type but totally negative 

CM . lepton asymmetry. 
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CM ' I. INTRODUCTION 

CM ' 

\ The success of big bang nucleosynthesis (BBN) is one of the most powerful pieces of evidence of standard big 



bang cosmology |Q. Roughly speaking, the predicted primordial abundances of light elements (D, 3 He, 4 He, and 7 Li) 
coincide with those inferred from observations for the baryon-to-photon ratio ri ~ 5 x 10 -10 . However, as observations 
are improved and their errors are reduced, a small discrepancy may appear Furthermore, rj is also determined by 
observations of small scale anisotropies of the cosmic microwave background radiation (CMB)[|[ [|, ||, 0], which may 
also cause the small discrepancy. Of course, the discordance may be completely removed as observations are further 
improved. However, it is also probable that such small discrepancies are genuine and suggest additional physics in 
BBN. 

These discrepancies are often eliminated if predicted primordial abundance of 4 He is decreased. Such a decrease is 
\ realized if there exists large and positive lepton asymmetry of electron type [Q . This is mainly because the excess of 
electron neutrinos shifts the chemical equilibrium between protons and neutrons toward protons, which reduces the 
O ,' predicted primordial abundance of 4 He. Note that this effect is much more effective than the corresponding speed 
up effect, that is, an increase of the Hubble expansion due to the presence of the chemical potential, which makes 
' the predicted primordial abundance of 4 He increase. However, large and positive lepton asymmetry of the electron 
^ , type is incompatible with small baryon asymmetry if we take account of the sphaleron effects, which convert lepton 
asymmetry to baryon asymmetry of the same order with the opposite sign [|| . This problem is evaded if one of the 
following three conditions is satisfied, that is, (a) lepton asymmetry is generated after electroweak phase transition 
but before BBN, (b) sphaleron processes do not work, (c) positive lepton asymmetry of the electron type is generated 
but no total lepton asymmetry is generated. 

The first condition was discussed in the context of neutrino oscillations || . In this case, large lepton asymmetry is 
generated through oscillations between active neutrinos and sterile neutrinos. In order to explore the second condition, 



one should note that the presence of large chemical potential prevents restoration of electroweak symmetry 10 . Based 
on this fact, large lepton asymmetry compatible with small baryon asymmetry was discussed pd| . The third condition 
is discussed by March- Russell et al. Jl^] . The Affleck-Dine mechanism produces positive lepton asymmetry of electron 
type but no total lepton asymmetry, that is, L e = —L^ > and L T = for some flat direction, which generates small 
baryon asymmetry due to thermal mass effects of sphaleron processes. 

In this paper we consider another possibility, which is something like the combination of (a) and (b). The Affleck- 
Dine (AD) mechanism produces positive lepton asymmetry of the electron type but totally negative lepton asymmetry 
by choosing an appropriate flat direction in the minimal supersymmetric standard model (MSSM) E3]. As an example, 
we identify the "ef L2L3" flat direction to be the AD field and consider the Affleck-Dine leptogenesis. Here subscripts 
represent the generations. Then, L e = —L^ = —L r = — Ltotal > is realized. The shift of the chemical equilibrium 
between neutrons and protons due to the positive chemical potential of electron neutrinos affects the results of BBN 
dominantly, while the speed-up effect caused by all the species of neutrinos is relatively negligible. After the Affleck- 
Dine leptogenesis. the AD field experiences spatial instabilities and deforms into nontopological solitons, Q-balls 
(L-balls) iBL Then, almost all the produced lepton numbers are absorbed into the L-balls |l6|, |l^|. If the 
lifetime of such L-balls is longer than the onset of electroweak phase transition but shorter than the epoch of BBN, 
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the large lepton asymmetry is protected from sphaleron effects and later released into the universe by the decay of the 
L-balls. On the other hand, small (negative) lepton numbers are evaporated from the L-balls due to thermal effects 
before the electroweak phase transition, which are transformed into small baryon asymmetry through the sphaleron 
effect. 

In our scenario we consider the Affleck-Dine mechanism and the subsequent Q-ball formation in the gauge-mediated 
supersymmetry (SUSY) breaking model. This is mainly because, in the gravity-mediated SUSY breaking model, the 
energy per unit charge of Q-balls is large enough to produce the lightest supersymmetric particles (LSPs) so that 
they will overdose the universe. Therefore, we do not consider the gravity-mediated SUSY breaking model. Since 
we assume that the AD field starts oscillating from the gravitational scale to produce large lepton asymmetry, the 
produced Q-balls are "new" Q or "delayed"- type [jlTj, depending on the sign of the coefficient of the one- loop 
correction to the effective potential. However, since the decay processes of the new type Q-balls are not completed 
before BBN, our scenario does not apply for them. Thus, we concentrate on delayed-type Q-balls in gauge-mediated 
SUSY breaking models. 

The rest of the paper is as follows. In Sec. II, we briefly review the Affleck-Dine mechanism and properties of 



Q-balls. In Sec. Ill, we discuss our mechanism to generate large lepton asymmetry compatible with small baryon 



asymmetry. Section [V is devoted to discussion and conclusions. 



II. AFFLECK-DINE MECHANISM AND Q-BALL FORMATION 

In this section we briefly review the Affleck-Dine mechanism and properties of Q-balls. In MSSM, there exist flat 
directions, along which there are no classical potentials in the supersymmetric limit. Since flat directions consist of 
squarks and/or sleptons, they carry baryon and/or lepton numbers, and can be identified as the Affleck-Dine (AD) 
field. In the following discussion, we adopt the "e c LL" direction as the AD field. In this case the AD field carries 
only the lepton number. 

These flat directions are lifted by supersymmetry (SUSY) breaking effects. In the gauge-mediated SUSY breaking 
model, the potential of a flat direction is parabolic at the origin, and almost flat beyond the messenger scale |l4|, [l?], [l9) , 



mj|*| 2 (|*| « Ms), 

"'-^Mjflogff) 2 (1*1 »M S ), W 



where is a soft breaking mass ~ 0(1 TeV), Mp is the SUSY breaking scale, and Ms is the messenger mass scale. 
Since gravity always exists, flat directions are also lifted by gravity- mediated SUSY breaking effects pQ], 



Vgrav — ^3/2 



l + iHog 



M 2 



l<*f, (2) 



where K is the numerical coefficient of the one-loop corrections and M is the gravitational scale (~ 2.4 x 10 18 GeV). 
This term can be dominant only at high energy scales because of small gravitino mass ^ 0(1 GeV). 

There is also the thermal effect on the potential, which appears at two-loop order as pointed out in Ref. . This 
effect comes from the fact that the running of the gauge coupling g(T) is modified by integrating out heavy particles 
which directly couple with the AD field. This contribution to the effective potential is given by 

V^-calT* log (3) 

where |c| ~ 1, and a w = g^/^ir represents the gauge coupling constant of the weak interaction since we consider e c LL 
direction. Though the sign of c depends on flat directions, it is irrelevant to our discussion since we assume that the 
zero-temperature potential dominates over the thermal effects. Note that a w should be replaced with a s for those 
flat directions which contain squarks. 

The lepton number is usually created just after the AD field starts coherent rotation in the potential, and its number 
density is estimated as 

n L {t osc ) ~ euj<j) 2 osc , (4) 

where e(^S 1) is the ellipticity parameter, which represents the strongness of the A term, and u> and cf> osc are the 
angular velocity and amplitude of the AD field at the beginning of the oscillation (rotation) in its effective potential. 
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Actually, however, the AD field experiences spatial instabilities during its coherent oscillation, and deforms into 
nontopological solitons called Q-balls (l4| [l5|, |l6) . When the zero-temperature potential V gauge dominates at the onset 
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El 110. 

of coherent oscillation of the AD field, the gauge-mediation type Q-balls are formed. Their mass Mq and size Rq are 
given by [El 



M Q ~ M F Q 3/4 , R Q ~ M^Q 1 / 4 . (5) 

From the numerical simulations |l7 [^the produced Q-balls absorb almost all the charges carried by the AD field 
and the typical charge is estimated as] 17fl 

with (3 fa 6 x 1CT 4 . 

There are also other cases where V grav dominates the potential at the onset of coherent oscillation of the AD field. If 
the coefficient of the one- loop correction K is negative, the gravity-mediation type Q-balls ( "new" type) are produced 
On the other hand, if K is positive, Q-balls do not form until the AD field leaves the V grav dominant region. 
Later it enters the V gauge dominant region and experiences instabilities so that the gauge-mediation type Q-balls are 
produced (delayed- type Q-balls) |L7[ . 

In our scenario described in the next section, the AD field starts to oscillate from the gravitational scale, i.e., 
<\>osc — M, which leads to the formation of new or "delayed" -type Q-balls. However, our scenario does not work for 
new type Q-balls because the produced Q-balls are large and do not decay before BBN. Hence, we concentrate on the 
delayed-type Q-balls below. Since the sign of K is in general indefinite and dependent on the model of the messenger 
sector in gauge-mediated SUSY breaking models, we assume that K is positive and delayed-type Q-balls are formed. 

When the AD field starts to oscillate in the V grav dominant region, where H osc ~ to ~ TO3/2, the lepton number 
is produced as ul — zmz/z&tsc- Since the delayed-type Q-balls are formed only after the AD field enters the V gauge 
dominant region for positive K 7 the charge of Q-ball is given by 



(7) 



with <f) eq ~ Mp/m^/2- Here the subscript "eq" denotes a value when the gauge- and the gravity- mediation potentials 
become equal. Thus the delayed-type Q-balls are formed at H eq ~ Mp/M . 

As we mentioned above, Q-balls absorb almost all the charges carried by the AD field. If we adopt the e c LL 
direction, all the lepton charges are confined in the Q-balls, namely, L-balls. Consequently, we must take out lepton 
charge from the L-balls through the evaporation, diffusion, and their decay. Part of the evaporated lepton charge is 
transformed into baryon charge by the sphaleron process, which accounts for the present baryon asymmetry. 

In the case of L-balls, they decay into leptons such as neutrinos via gaugino exchanges. The decay rate of Q-balls 
is bounded as [E3| 



dQ 



dt 



<^ (8) 
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where A is a surface area of the Q-ball. For L-balls, the decay rate is estimated as a value of the order of the upper 
limit. 

According to Refs.jT?], 0, we evaluate the evaporation rate of L-balls, which is given by (24) 

Covap = = — ^(MQ — Aiplasma)? l2 47ri?Q, 

~ -AtTK^1qT 2 R q for flQ > ^plasma, (9) 

where /iq and /i p i aS ma are chemical potentials of the Q-ball and plasma, and the coefficient n 1 includes statistical 
and other numerical factors. The chemical potential of the Q-ball is given as /iq ~ u since the energy of the cf> particle 
inside the Q-ball is u. At m^, large numbers of the scalar particles building up Q-balls are in the plasma, which 
implies k ~ 1. On the other hand, at m$, the evaporation from Q-balls is suppressed by the Boltzmann factor. 
In the case of L-balls, the main process of the evaporation is 4>4> — * 11 through wino or bino exchange, which yields 
k ~ a^T 2 1 'mi at m^. 
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However, if the charge transport is not effective enough, the evaporated lepton charges in the "atmosphere" of 
the L-ball will establish chemical equilibrium there. In this case, the dissipation of the charge is determined by the 
diffusion. The diffusion rate is estimated as p5| , 

Cdiff = = -^DRq^qT 2 

~ ~4nDT 2 , (10) 
where the diffusion constant D of relativistic sleptons and leptons in a hot plasma is given by D ~ a/T with 



In short, the time scale of the charge transportation is determined by the evaporation rate when 
ICcvapl <"fCdiffI7and by the diffusion rate when |Cova P | > | Cdiff I - 

The amount of the evaporated charges can be estimated by integrating Eqs. (^) and ([l0|) in the course of the 
evolution of the universe. When the AD field starts to oscillate at the gravitational scale, its oscillation energy is 
comparable to the total energy of the universe. Therefore, the energy of the universe will be dominated by the AD 
condensate or Q-balls soon after the reheating and the universe continues to be matter dominated. The thermal history 
of the universe is rather involved because radiation comes from both decays of an inflaton and Q-balls. However, in 
fact, we have only to consider two cases where the cosmic temperature decreases monotonically. 



III. LARGE LEPTON ASYMMETRY FROM L-BALL 



In this section we give a detailed explanation of our scenario. Our goal is to generate small baryon asymmetry and 
positive large lepton asymmetry of the electron type simultaneously. In general, however, this is difficult to accomplish 
because the chemical equilibrium induced by the sphaleron transition forces the baryon and the lepton asymmetries 
to be of the same order with opposite sign ||. Hence we must get over two problems : (i) how to protect large 
lepton asymmetry from being converted to baryon asymmetry by the sphaleron process, and (ii) how to reconcile the 
opposite sign of baryon and lepton asymmetries. 

We show that these two obstacles can be evaded by considering the Affleck-Dine lcptogenesis and the subsequent 
L-ball formation using e c LL direction. First of all, we give the outline of our scenario and the solution to the problem 
(i). Large lepton asymmetry can be generated if the A terms, which make the AD field rotate in the effective potential, 
originate from some Kahler potential with vanishing superpotential. Then the AD field starts to oscillate with large 
initial amplitude 4> osc ~ M and ellipticity e ~ 1. As spatial instabilities grow, delayed- type L-balls are formed 
and absorb almost all charges carried by the AD field. It is essential to our scenario that the lepton asymmetry 
confined in the L-balls is kept from the sphaleron process. However, a small part of lepton charges confined in the 
L-balls are evaporated due to thermal effects. Thus, lepton charges evaporated until the electroweak phase transition 
(T^,Tc ~ 300 GeV) AQ ew are partly converted to baryon asymmetry through the sphaleron process, which explains 
the present small baryon asymmetry. On the other hand, large lepton asymmetry comes out through the decay of the 
remnant L-balls after the electroweak phase transition, which must be completed before BBN. Thus the small ratio 
AQew I Q is the source of hierarchy between the baryon and the lepton asymmetries. 

Next we give a solution to the problem (ii), that is, the sign of the lepton asymmetry. What we want to generate is 
positive baryon asymmetry and positive lepton asymmetry of electron type. However, the sphaleron process converts 
positive lepton asymmetry into negative baryon asymmetry. To surmount this problem, we adopt the e\LiL^ direction 
as the AD field, which leads L e = —L^ = —L r = — Ltotai > 0. Thus the positive lepton asymmetry of the electron 
type is generated, whilst total lepton asymmetry is necessarily negative in order to have positive baryon asymmetry 
through the sphaleron transition. At the epoch of BBN, charged leptons except electrons have already disappeared 
through decay and annihilation processes. Also, because of the charge neutrality of the universe, lepton asymmetry 
stored in electrons are comparable to baryon asymmetry, which is rather small. Thus, there can exist large lepton 
asymmetry only in the neutrino sector. For later use, we define the degeneracy parameter as the ratio of the 
chemical potential to the neutrino temperature. The presence of chemical potentials speed up the universe, which 
leads to an increase in the n/p ratio. However, its effect is negligible in comparison with the effect of the shift of 
chemical equilibrium between protons and neutrons due to the chemical potential of electron neutrino in the case of 

\U = \^\ = \U\- 

Now we give a quantitative estimate for our scenario. We assume that the zero-temperature potential dominates, 

(2) 

i.e., Vg aU g e ^> Vj, , at the formation of the delayed- type L-balls with H ~ H eq : 

a 2 w T* q < M*, (11) 

where T eq is the temperature of the universe just before the delayed-type L-balls are formed. As shown below, this 
constraint is automatically satisfied for the cases we consider. 
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The delayed-type L-balls must decay before BBN, 



dQ 



dt 



< lsec, 



which leads to the constraint 



m 3/2 
lOMcV 



> 



Mp \ 



4/5 



lOTeV 



(12) 



(13) 



Here Eq. (g) is used. In order to estimate the baryon and the lepton to entropy ratio, it is necessary to evaluate the 
entropy production by the decay of the L-balls. The decay temperature of the L-balls, Td, is given by 



T d = 



90 



1/4 



M- 



48tt 



~ 1.3MeV 



/ M F \ / TO3/2 \ 5 / 2 



VlOTcvJ V 



lOMeV 



) 



(14) 



where g* = 10.75 counts the total number of effectively massless degrees of freedom. 

Now we turn to an account of the total evaporated charge, AQ, and the evaporated charge at temperatures above 
the electroweak phase transition, AQ ew . In fact we have only to consider the following two cases. In the other cases, 
the temperature during the presence of L-balls does not exceed Tq so that the evaporated lepton numbers are not 
converted into baryon numbers. 

First we consider the case that the delayed-type L-balls are formed before the reheating and decay after that (case 
A). This is realized if the following two conditions are satisfied, 



M F > T RH , 



(15) 
(16) 



Then the temperature at the L-ball formation is given by T eq = \JMpTr R , which automatically satisfies the require- 
ment (|TT|) . The temperature of the universe is approximately given as 



1/3 



' (T 2 RH MH) 1 for T RH < T, 
(TjMH) 1 " 

, Vhm 



for T p <T < T RH , 
for T d <T< T p , 
for T < T d , 



(17) 



where T p = (T R hT£) denotes the temperature when the radiation derived from the decay of the L-balls dominate 
over those derived from the inflaton. Note that the cosmic temperature decreases monotonically in this case. The 
condition that the chemical equilibrium induced by the sphaleron transition are well established is given by 



T eq = \J MpT RH > T c . 

With the use of Eq. (0) , the evaporation rate with respect to the temperature is estimated as 



(18) 



dQ 
dT 



4tt 



M F T 3 



,1/4 



4-7ra„ 



47T- 



T 2 M 

1 RH IU 



mAMpT 



Q 



1/4 



evap 



MpVTTrh 
MT 3 / 2 



Q 



1/4 



for 7710, T R H < T < T eq , 

for T RH <T <m (j „Teq, 
for m^, T p <T < T RH , 



4.5^Qi/ 4 



— Q 1 / 4 for T v < T < m#1 T RH , 



(19) 



M F T 3 



Ana 



2 T 2 M 
'm 2 M F T 



)l/4 



for m^, T d <T < T p , 
for T d <T < TTifp , T p . 
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On the other hand, the diffusion rate with respect to the temperature is given by 

T 2 RH M 



dT ) 



Ana- 



Ana 



diff 



Ana- 



M 

T l/2 T 3/2 

TjM 



r 4 



for T RH <T <T eq , 
for T p <T< T RH , 

for T d <T< T p . 



(20) 



By integrating Eqs. and (£0|), the evaporated charges AQ and AQ ew are found to be the same order, and given 
by 



AQ ~ AQ e 



AnaT RH M 



3 

87ra- 



M 



Ana TjM 

m? A 



for T RH < m,p < T eq , 
for Tp < 7770 < Trh, 
for T d < < T p , 



(21) 



where we have used Tc ~ m^. 

Next we consider the case where delayed-type L-balls are formed after the reheating and the temperature decreases 
monotonically (case B). This is realized if the following conditions are satisfied, 



T RH > M F , 
M F > (T RH T d/ 



3\ 5 



(22) 
(23) 



1 /3 

Then the temperature at the Q-ball formation is given by T eq — \M f /Trh) , which again satisfies the requirement 
(|ll|). Though the time evolution of the cosmic temperature is the same as case A, the requirement for the sphaleron 
process to work now reads 



T, 



(Mp/T, 



RH 



,1/3 



>T C . 



The evaporated charges AQ and AQ ew can be estimated similarly and given by 

M 



8na- 



for T p < 77i0 < T e , 



AQ ~ AQ e 



J tti^Trr 

forT (i <m <r p 



(24) 



(25) 



Finally we estimate the baryon (lepton) to entropy ratio, using the results derived above. The baryon to entropy 
ratio is then given by 



n B 
s 



m 3/2 M' fog.T* AQ e 



2 T d AQ e 
23 



3/2 



m-3/2 



4 x 10" 



ITeV 



/m 3/a w / Trh V ( M F \ ' 6 

VIGeV/ ^lOGeVy V10 6 GeVy 

for Tflij < 7770 < T e(? (caseA) 

3 x 10 6 GeV / 

for T p < 7770 < T eq (caseB) 



dXW \TTcVJ VlGcvJ llO^GevJ 



where we have used Eqs. (@), @), ©, and (||). 
way, the lepton number to entropy ratio is given by 



(26) 



Also we have assumed the maximal CP violation. In the same 



s 



T, 



Am 



3/2 



7 



-2 



which yields |7|] 

£„ e ~ -10 x — 



The allowed regions for m 3 / 2 7 Mp, and is shown in Fig. |], where the baryon to entropy ratio takes the value 
required from BBN, 

10- 11 <!^<10- 10 . (29) 
~ s ~ 

Here we adopt a rather loose constraint because of the uncertain CP phase. As can be seen from Fig. [l], there are 
two allowed regions : (i) m 3/2 ~ 0.1 - lGeV, M F ~ 10 5 - 10 6 GeV, and T RH ~ 1 - 10 3 GeV, (ii) m 3/2 - 0.1 - lGeV, 
Mp ~ 10 6 GeV, and Trh ~ 10 6 — 10 9 GeV. Roughly speaking, the regions (i) and (ii) correspond to the cases A and 
B respectively. 

Also we plot the contours of the degeneracy of electron neutrinos in Fig. [^, which shows that the large and positive 
lepton asymmetry of electron type can be generated in our scenario. For reference, the present constraint of £„ e by 
the analyses of BBN and CMB data is given by |2J| , 

-0.01< £„ e <0.22. (30) 

Thus, our scenario can generate both small baryon asymmetry and positive large lepton asymmetry of electron type 
at the same time by virtue of the AD leptogenesis and subsequently formed L-balls. 

IV. DISCUSSION AND CONCLUSIONS 

In this paper we have proposed a scenario which accommodates small baryon asymmetry and large lepton asymmetry 
simultaneously. The large lepton asymmetry is generated through the Affleck-Dine mechanism and almost all the 
produced lepton charges are absorbed into L-balls which are formed subsequently. Thus, most of the produced lepton 
numbers do not suffer from the sphaleron process. Only a small fraction evaporated from the L-balls due to thermal 
effects is converted into baryon asymmetry, which is responsible for the present baryon asymmetry. 

As a concrete example, we consider positive large lepton asymmetry of the electron type. The excess of electron 
neutrinos shifts the chemical equilibrium between protons and neutrons toward protons so that the predicted primor- 
dial abundance of 4 He is decreased, which often gives a solution to the discrepancy of BBN itself or that between 
BBN and CMB. However, the sphaleron process converts lepton asymmetry into baryon asymmetry with the opposite 
sign. To circumvent this problem, we identify the e\L2L 3 flat direction to be the AD field. Then the Affleck-Dine 
leptogenesis can generate positive lepton asymmetry of the electron type but totally negative lepton asymmetry, which 
is converted into positive baryon asymmetry. Of course, one should notice that by use of another flat direction such 
as e^LiLs, we can obtain negative lepton asymmetry of the electron type and also total negative lepton asymmetry. 

Recently, it was pointed out that complete or partial equilibrium between all active neutrinos may be accomplished 
through neutrino oscillations in the presence of neutrino chemical potentials, depending on neutrino oscillation pa- 
rameters fjof . In the case of partial equilibrium, our scenario needs no change. Only complete equilibrium can spoil 
our scenario. Even if neutrino oscillation parameters lead to complete equilibrium, our scenario may still work since 
it is possible that the L-balls decay just before BBN and the complete equilibration cannot be attained, which needs 
further investigation. 
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FIG. 1: The allowed region for m 3 / 2 , Mf and Trh, where our scenario succeeds and the baryon to entropy ratio satisfies 
the following bounds : 10 -11 ^Sna/s^S 10~ 10 . Note that there does not exist any upper bound on Trh from the gravitino 
problem (li], ^] , since the L-balls dominate the universe and their decay temperature is rather low. The two separate allowed 
region roughly corresponds to the cases A and B discussed in the text. 
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FIG. 2: The contours of the electron neutrino degeneracy are shown. The trapeziform area between two solid lines represents 
the allowed region where our scenario works and the baryon to entropy ratio satisfies the bounds: 1CP 11 ^ ub/s S 1CF 10 . The 
contours represents £„ e = 0.005, 0.01, 0.02, 0.06, 0.1 from top to bottom. 



